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Abstract 



^ ' The fundamental properties of biquaternions (complexified quaternions) are presented including 

several different representations, some of them new, and definitions of fundamental operations such 
as the scalar and vector parts, conjugates, semi-norms, polar forms, and inner and outer products. 
The notation is consistent throughout, even between representations, providing a clear account of 
the many ways in which the component parts of a biquaternion may be manipulated algebraically. 
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1 Introduction 



It is typical of quaternion formulae that, though they be difficult to find, once found they are 
immediately verifiable. 

J. L. Syngc (1972) [43, p34] 

The quaternions are relatively well-known but the quaternions with complex components (complexified 
quaternions, or biquaternions 1 ) are less so. This paper aims to set out the fundamental definitions 
of biquaternions and some elementary results, which, although elementary, are often not trivial. The 
emphasis in this paper is on the biquaternions as an applied algebra - that is, a tool for the manipulation 
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1 Biquaternions was a word coined by Hamilton himself [21] and [17, § 669, p664]. (The word was used 18 years later by 
Clifford [5] for a different concept, which is unfortunate.) 
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of algebraic expressions and formulae to allow deep insights into scientific and engineering problems. It is 
not a study of the abstract properties of the biquaternion algebra, nor its relations with other algebras. 
Throughout the paper 'quaternion' means a quaternion with real components (a quaternion over the 
reals, R), and 'biquaternion' means a quaternion with complex components (a quaternion over the field 
of complex numbers, C). We denote the set of quaternions by H, and the set of biquaternions by B. H 
is, of course, a subset of B. 

Some of the material in the paper is based on the book by Ward [45, Chapter 3] which is one of 
the few sources of detail on the biquaternions. We have not followed Ward's notation in this paper, 
preferring instead a scheme based on bold or plain symbols without hats and underlines. 

We have also drawn upon the paper by Sangwine and Alfsmann [35, 36] which sets out comprehen- 
sive results on the divisors of zero, and their subsets the idempotents and nilpotents. Sangwine and 
Alfsmann's paper uses the same notations as this paper (having benefited from access to this paper in 
draft). 

The quaternions themselves (with real elements) are well-covered in various books, for example [1, 2, 
27, 28, 29]. Hamilton's works on quaternions were published in book form in [17, 18, 14], and many are 
also now available freely on the Internet in various digital repositories. The paper by Coxeter [6] is also 
a useful source. 

We begin by setting out various ways in which quaternions and biquaternions may be represented. We 
start with representations for quaternions, for reference, and because these representations are generalized 
for biquaternions. In what follows, we use notation as consistently as we can. In particular, the complex 
operator usually denoted by i (or in electrical engineering j) is represented in this paper as / in every 
case. This is to keep the complex operator distinct from the first of the three quaternion operators i, 
since it is independent. The independence of i and / is perhaps the most fundamental axiomatic aspect 
of the biquaternions that must be understood. Bold symbols denote vectors and bivectors 2 , whereas 
normal weight symbols denote scalar or complex quantities or quaternions. 

Throughout the paper 3 we use the term norm or the more specialized term semi-norm, both denoted 
to mean the sum of the squares of the components of a quaternion or biquaternion, and modulus, 
denoted to mean the square root of the norm and thus the Euclidean magnitude. This is not 
universally accepted terminology, many sources using norm where we use modulus. However, our usage 
is consistent with several authors who have written on quaternions, including Synge [43] and Ward [45], 
but it does require care when using statements made about norms in other sources. 

Many of the concepts given in this paper are implemented in numerical form in a MATLAB® toolbox 
[37] which two of the authors first developed in 2005 and are built upon in a toolbox for handling linear 
quaternion systems, first developed in 2007 [10]. The toolbox [37] was essential to the development of 
this paper and the results presented within it, since otherwise, errors in algebra would go unnoticed. In 
many cases we have established results first by using the toolbox, and then derived the algebraic proofs 
or statements which appear here. 

2 Quaternions 

Classically, quaternions are represented in the form of hypercomplex numbers with three imaginary 
components. In Cartesian form this is: 

q = w + xi + yj + zk (1) 

where i, j and k are mutually perpendicular unit bivectors obeying the famous multiplication rules: 
i 2 = j 2 = k 2 = ijk = — 1, discovered by Hamilton in 1843 [19], and w, x, y, z, are real. Quaternions 
are generalized to biquaternions by permitting w, x, y and z to be complex, as discussed in § 3, but in 
this paper we reserve the symbols w, x, y, z for the real case. A quaternion with w = is known as a 
pure quaternion (Hamilton's terminology, but still used and widely understood). 

2 Bivcctors represent directed areas, and are explained in Table 3 and §3.6. 
3 Except in equation 11, where we use the notation of the cited reference. 

4 Classic texts often refer to the operators i, j and k as vectors, a misconception that has caused considerable confusion 
over many years, but is understandable, since it could not be cleared up without the concept of geometric algebra and 
bivectors. We discuss this in § 3.6. 
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The conjugate of a quaternion is given by negating the three imaginary components: q = w — xi — 
yj — zk. It is easily shown that qp — pq for general quaternions p and q. Indeed the formula may 
be generalized to more than two quaternions (the generalized formula was first noted by Hamilton [20, 
§20, p238], and was also included in [28, p60, §31]): pqrst = tsfqp. The quaternion conjugate may be 
expressed in terms of multiplications and additions [11, Theorem 11] using any system of three mutually 
orthogonal unit pure quaternions (here i,j,k): 

q = -^(q + iqi+ jqj + kqk) 

Similar formulae, based on involutions [9, 11], exist for extracting the four Cartesian components of a 
quaternion 5 [41]: 

(q - iqi - jqj - kqk) 



1 

4 
1 

Ti 
l 

ij 
l 

4k 



(q - iqi + jqj + kqk) 

(2) 

(q + iqi - jqj + kqk) 
(q + iqi + jqj - kqk) 



The norm of a quaternion is given by the sum of the squares of its components: ||g|| = w 2 + x 2 + y 2 + 
z 2 , \\q\\ G K. It can also be obtained by multiplying the quaternion by its conjugate, in either order since 
a quaternion and its conjugate commute: ||g|| = qq = qq. The modulus of a quaternion is the square 
root of its norm: |g| = ylMI- 

Every non-zero quaternion has a multiplicative inverse 6 given by its conjugate divided by its norm: 

The quaternion algebra H is a normed division algebra, meaning that for any two quaternions p and 
q, \\pq\\ = \\p\\ \\q\\, and the norm of every non-zero quaternion is non-zero (and positive) and therefore 
the multiplicative inverse exists for any non-zero quaternion. 

Of course, as is well known, multiplication of quaternions is not commutative, so that in general for any 
two quaternions p and q, pq ^ qp. This can have subtle ramifications, for example: {pq) 2 = pqpq ^ p 2 q 2 ■ 

Alternative representations for quaternions are given in Table 1, expressed in terms of the Cartesian 
form given above in equation 1, and in selected cases, in terms of other representations given here, fj, is 
a unit pure quaternion and is known as the axis of the quaternion. It expresses the direction in 3-space 
of the vector part 7 , V(g). Hamilton himself showed that any unit pure quaternion is a square root of 
— 1 [22, pp203, 209][17, § 167, pl79] and a proof is also given in [9, Lemma 1]. This is why we call the 
form a + fib the 'complex' form, since it is isomorphic to a complex number a + lb. This means that 
the modulus and argument of the quaternion are identical to those of a + lb (one could think of them 
as having similar Argand diagrams, where in the case of quaternions, the Argand diagram represents 
a plane section of 4-space defined by the 'axis' \x and the scalar quaternion axis - along which w is 
measured) . 

The polar form of a quaternion is analogous to the polar form of a complex number, with one 
exception. The argument, 9, is confined to the interval [0,7r) because the modulus of the vector part is 
always taken to be positive (there is no convenient way to define an orientation in 3-space which would 
permit the sign of the vector part to be determined) . If we negate the argument of the exponential in the 
polar form, therefore, the negation is conventionally applied to the axis, \x and not to the argument 9. 
When 9 is computed numerically, the result is always in [0, n) because we have to use the (non-negative) 
modulus of the vector part to compute it (using an atan2 function, typically). 

The Caylcy-Dickson polar form [38, 39] has a complex modulus and a complex argument (both in 
fact are degenerate quaternions of the form w + ix 7 isomorphic to complex numbers). 



5 These formulae, and that for the quaternion conjugate, generalize to biquatcrnions, even with mutually orthogonal 
unit pure biquatcrnions in place of i, j and k, although this latter point has not been thoroughly checked. 
6 This does not apply to biquatcrnions. 

7 We use the term 'vector part' throughout this paper to mean that part of the quaternion consisting of the three 
components containing i, j and k. It does not necessarily correspond to the concept of a 3-space vector, since it could be a 
vector, bivector or a combination of both, using the language of geometric algebra described later. The term 'vector part' 
is well-established in the literature and, lacking a good alternative which would be readily understood, we retain it. 
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Designation 


Representation 


Details (w, x, y, z, a, b, r, 9 £ IR, q £ H, 
fi G H,S(/x) = 0,/x 2 = -1) 


Cartesian 


q = w + xi + yj + zk 




Scalar + vector 


q = S(q)+V(q) 


S{q) = W 

v \q) — xi -\- yj + zk 


'Complex' form 


q = a + lib 


a = w 

b= V^ 2 + V 2 + z 2 = |V(g)| 
M — v \1) 1 \ v \Q)\ 


Cayley-Dickson 


q = (w + xi) + (y + zi)j 


This multiplies out to the Cartesian 
representation. 


Polar form 


q = r cxp (fiff) = r (cos 9 + fj, sin 9) 


r = \q\ 
r cos 9 = a 
r sin 9 = b 


Cayley-Dickson 


q = Acxp(Bj) 


A,B eU with y = z = 


polar form 




(isomorphic to C) 

See [38, 39] for formulae defining A and B. 



Table 1: Representations for quaternions. 



3 Biquaternions 

To generalize the quaternions to biquaternions we simply permit the four elements to be complex rather 
than real, thus giving us the Cartesian representation: 

q = W + Xi + Yj + Zk (3) 

where i, j and k are exactly as in § 2 and W, X, Y, Z, are complex. This generalization was first studied 
by Hamilton himself [2f] and [17, §669, p664], and was also discussed by Cay ley [44]. 

In equation 3 each of the four elements is of the form W = %t(W) + I?s(W), where I 2 = — 1 is the 
usual complex operator, distinct from i. Axiomatically, / commutes with the three quaternion operators 
i, j and fc, that is il = li. jl = Ij and kl = Ik. Since reals commute with the three quaternion 
operators, so do all complex numbers. It is also worth noting that any quaternion with only one non- 
zero component in the vector part, for example w + ix is not a complex number, but a degenerate 
quaternion. Such numbers are isomorphic to the complex numbers, but it is best not to confuse them 
with complex numbers when working with biquaternions. 

Some familiar rules of algebra apply to biquaternions. Since the real and imaginary parts are quite 
separate from the concept of the four quaternion components, real and imaginary parts may be equated 
just as when working with complex equations. However, some of the elementary properties of quaternions 
become non-elementary when the quaternions are complexified. For example, generalizations of the norm 
and modulus of a quaternion are complex in general for biquaternions, and so is the argument in one of 
the polar forms. We discuss each of these non-trivial properties in a later section. 

The existence of complex generalizations of the norm, modulus and inner product, yields a problem 
of terminology, since conventionally these quantities are real, and satisfy properties that their complex 
generalizations cannot (the triangle inequality \\p + q\\ < \\p\\ + \\q\\, for example, requires ordering, but 
complex numbers lack ordering, hence the triangle inequality cannot be applicable to a 'norm' with a 
complex value). Rather than invent new terms, we use the existing accepted terms (with the exception 
of norm, where we substitute semi-norm, for reasons discussed in §3.3), but caution the reader that 
because these quantities are complex, they cannot be assumed to satisfy all the usual properties of their 
conventional real equivalents. In taking this approach, we are following Synge [43, p9], and to some 
extent Ward [45], both of whom use the term norm to refer to a complex generalization. Synge also 
refers to a scalar product with a complex value, and Ward uses the term inner product for a complex 
generalization of the concept in his section on the Minkowski metric [45, § 3.3, p 115]. 

Although a biquatcrnion commutes with its quaternion conjugate, and complex numbers commute 
(including with their complex conjugates), somewhat surprisingly, a biquaternion does not necessarily 
commute with its complex conjugate (Proposition 1 in §3.1). 
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Alternative representations for biquaternions are shown in Table 2, expressed in terms of the Cartesian 
form given above in equation 3, and in selected cases, in terms of other representations given here. In 



Designation 


Representation 


Details: w a , x a , y a , z a g R, a £ {r, £}, 
W,X,Y,Z,A,B,R,Q G C, 
q r ,qi,Q,^ ettgel, 

£ d M £ 2 — 1 \"XA 


Cartesian 


q = W + Xi + Yj + Zk 




Scalar + vector 


q = S(q) + V(q) 


o i y i — i' v 

V \q) =a« + i J + Zk 






A = W = 8(q) 


'Complex' form I 


q = A + £B 


£B = V(q) 

B = \/ X z + Y 2 + Z 2 = \\{q)\ 
£= (Xi + Yj + Zk)/B 


'Complex' form II 


q = q r + Iq, = Rfa) + m(q) 


q r = w r + x r i + y r j + z r fc 

n ■ 1 1 ) - 1 T* ■ t 1 T J - 1 1 . Jt> 

Hi — IDiTXitT i/jj l *»*«■ 






w r = 3?(W^),^ = 


Expanded form 


q = w r + x r i + y r j + z r k + 
I(wi + Xi% + yd + Zik) 


x r = ^(X),^ = 3(X) 

wr = »(r),»i = 9f(r) 

z r = 5R(Z),Zj = $t(Z) 


Hamilton polar form 


q = i?cxp(£9) 
= R (cos 6 + £ sin 9) 


R=\q\ 

A = RcosS 

B = RsinQ 






* = tan" 1 (q- 1 q t ) 


Complex polar form 


q = Qexp(I*) 
= Q (cos* + / sin*) 


Q = q/ exp (/*) = q exp (-/*) 
q r = Q cos * 
(ft = Q sin * 



Table 2: Representations for biquaternions. 



the scalar/vector form, the scalar part is complex, and the vector part is a pure biquaternion. 

Complex form I corresponds to the 'complex' form in Table 1. The differences are that A and B are 
now complex, whereas a and 6 were real, and the imaginary unit £ is a pure biquaternion root of —1 of 
the form £ = bfi + dlv, where b 2 — d 2 = 1 and fi and v are mutually perpendicular unit pure quaternions 
(themselves also roots of —1) [33, 34]. Note that B can vanish - as discussed in § 3.5.2 and in detail in 
[35, 36, § 2] - although its components (X, Y and Z), and therefore q, do not vanish. £ can be thought 
of as a complex axis, in the sense that it has real and imaginary parts which each define directions in 
3-space. The geometric interpretation of biquaternions is discussed further in §3.6. 

Complex form II has quaternions in the real and imaginary parts, and is perhaps the most obvious 
representation for a biquaternion other than the Cartesian form. It is related to the complex polar form 
described below. 

In the expanded form, the biquaternion is represented as a complex number with quaternion real and 
imaginary parts expressed in Cartesian form. 

The polar form of a quaternion depends on Euler's formula cxp(/(9) = cos 8 + 1 sin 8 which generalizes 
by replacing the complex operator / with any root of —1. In the case of quaternions, the set of unit pure 
quaternions provides an infinite number of roots of —1 and the general polar form rexp(/x#), as given 
in Table 1 is therefore a straightforward extension of Euler's formula. In the case of biquaternions there 
are two possibilities for the root of —1: the complex root /, or any one of the biquaternion roots of —1 
defined in § 3.4. Thus we have two possible fundamental polar forms. 

The first ('Hamilton') polar form generalizes the polar form in the real case. R, the 'modulus' of the 
biquaternion, is complex; £ is a biquaternion root of — 1; and 0, the 'angle' in the exponential, is also 
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complex 8 . The interpretation of complex angles is discussed further in §3.2. 

In the second ('complex') polar form, the standard complex operator I serves as the root of —1 in the 
exponential, but in this form, the 'argument' of the exponential, *, is a quaternion, and the exponential 
is scaled by a quaternion 'modulus' Q. We thus have a polar form with a 'modulus' and 'argument' 
in H. Note that it is not possible to find Q by the obvious direct route of Q = \J q% + qf because 
q r = Qcos* and squaring this gives Qcos* Qcos* and not Q 2 cos 2 *. However, it is the case that 
cos 2 * + sin 2 * = 1, as would be expected. Further, note that care is needed in computing the inverse 
tangent in order to find *: it is important that the real part is divided on the left. This is because 

tan* = (Qcos*) -1 (<2sin*) = cos -1 * Q -1 <3 sin* = cos -1 * sin* 
whereas with a division on the right we have: 

tan* ^ (Qsin*)(Qcos*) -1 = Q sin* cos -1 * Q -1 

and Q and its inverse are at opposite ends of the product and do not cancel. (Note that cos* and 
sin* commute, so their quotient is the same whether divided on the left or right.) The exponential is a 
biquaternion because of the presence of I and it has a complex modulus. The modulus of this complex 
modulus is 1. We discuss both polar forms further in §3.1 in the context of conjugation. Finally, note 
that in the complex polar form, Q and the exponential do not commute. It is therefore possible to define 
a variant by placing Q on the right of the exponential. The variant is related to the form in Table 2 by 
the conjugate rule. 

De Leo and Rodrigues [7, 8] discussed polar forms of biquaternions but apparently did not see that 
there were two simple polar forms as here, with different imaginary units. Instead they described a single 
polar form containing the product of two exponentials. We can do this with either of our polar forms, 
representing the 'modulus' in each case in its own polar form. Thus the 'Hamilton' polar form can be 
written: 

q = R cxp (£9) = r cxp (70) exp (£0) (4) 

where r, 4> £ R are the modulus and argument of i?, the complex 'modulus' of q. Note that, because I 
and £ commute, the two exponentials commute, and it is possible to write this polar form as: 

q = r exp(/</> + £6) 

where the argument of the exponential is now a biquaternion with I<j) as scalar part and £9 as vector 
part. (In general, with quaternions as well as biquaternions, e p e q ^ e pq because of non-commutativity.) 

Similarly, the 'complex' polar form can also be written in this way, expanding the quaternion 'mod- 
ulus', Q, into the standard polar form of a quaternion as given in Table 1: 

q = Q cxp (/*) = r exp (/Lt#) exp (/*) (5) 

where fj. £ M and r,6 £ R. Notice that the single real modulus r is the same in each case, but the various 
'angles' are different in value and type (0, 8 £ R, 9 £ C, * £ H). The real modulus, r, is the absolute 
value of the square root of the semi-norm as discussed in § 3.3. In this case, the arguments of the two 
exponentials (specifically \x and *) do not commute, hence the two exponentials cannot be combined by 
adding fi9 to /*, and neither can the order of the exponentials be changed. 

We can usefully combine the 'complex' representation of a quaternion from Table 1 with 'complex' 
form II from Table 2 to give the following representation, which was used in [34, 33] to derive the 
biquaternion roots of —1: 

q = (a + p0) + 7(7 + u8) (6) 

in which \x and v are real pure unit quaternions, and a, /3, 7 and S are real. The four real coefficients 
in this representation may be related to the coefficients in the other representations given above. For 

8 The cosine and sine of a complex angle are simply defined in terms of Euler's formula as cos 2 = i (e Iz + e~ Iz ) and 
sin z = — (e Iz — e~ IzS j ; or, if we write z = x+Iy: cos z = cos x cosh y — I sin x sinh y and sin z = sin x cosh y+I cos x sinh y 
[4, See complexes]. 
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example: a + 17 = A = W. The correspondence between £, and /x and u is not so simple. Equating the 
vector part of equation 6 with the vector part of 'complex' form I, we get: 

HP + IuS = £B = £y/X 2 + Y 2 + Z 2 

and we can see that dividing + IuS by its (complex) modulus B, will yield £ provided that B does 
not vanish, as discussed in § 3.5.2. 

We may relate each of the terms in equation 6 to a concept from geometric algebra [24, 40], as shown 
in Table 3 with equivalent definitions based on various representations from Table 2. 



Geometric algebra concept 


Term in equation 6 


Terms from Table 2 


Scalar - undirected quantity 


a 


ft(W)= n(A)= S(q r ) = w r 


Bivector - directed area 




R(££)= »(V(g))= V(«r) 


Vector - directed magnitude 


IuS 


J9(€S)=J9f(V(?)) = JV( ?i ) 


Pseudoscalar - undirected volume 


II 


ns(W)= ns(A)= is( qi )=iwi 



Table 3: Correspondence between geometric algebra concepts and biquaternion components. 

These equivalences are given by Ward [45, § 3.2, p 112] and they are discussed in more detail in § 3.6. 
Note carefully that geometric vectors are represented by imaginary pure quaternions, and that real pure 
quaternions are bivectors. This is because the product of two perpendicular vectors must yield a bivector. 
The product of two bivectors gives a bivector (and a scalar, unless the bivectors are perpendicular). 

The representation in equation 6 was used in [34, 33] to derive the solutions of the equation q 2 = — 1 
(that is the biquaternion roots of —1) when q is a biquaternion, and it was shown that the solutions 
required a = 7 = 0, /x _L u, and /3 2 — S 2 = 1. Clearly any of the other representations given in this 
section would not permit this result to be expressed so clearly, nor to be easily derived. We return to 
this result in Theorem 2. 

The biquaternion algebra B is not a division algebra because non-zero biquatcrnions exist that lack 
a multiplicative inverse. The set of such biquaternions is known as the divisors of zero [35, 36] and is 
defined in § 3.5. 

3.1 Conjugates 

The conjugate of a biquaternion may be defined exactly as for a quaternion by negating the vector part. 
Thus we have 

q = W-Xi-Yj-Zk = S(q)-V(q)=q^ + Iq- l (7) 

We call this the 'Hamiltonian' or quaternion conjugate, in agreement with Synge [43, Equation 3.4, p8]. 
The conjugate rule qp = pq and its generalization to more than two quaternions applies equally to 
biquatcrnions. Similarly, a biquaternion commutes with its quaternion conjugate, and the product of 
the two is the semi-norm, as discussed in § 3.3. 

However, there is another possible conjugate - that obtained by taking the complex conjugate of the 
complex components of the quaternion. We denote this complex conjugate by a superscript star, which 
is a common convention with complex numbers. Thus the complex conjugate is given by: 

q* = W* + X*i + Y*j + Z*k = S(q)* + V(q)* = q r - Iq, (8) 

Our definition in equation 8 agrees with that of Synge [43, p8, equation 3.4] and appears to be the obvious 
way to define the complex conjugate, but it differs from that of Ward [45] who defines a complex conjugate 
which is a quaternion conjugate with complex conjugate components. There is no way to define the 
complex conjugate in terms of additions and multiplications as can be done with the quaternion conjugate 
— if a means existed, it would also apply to complex numbers because a degenerate biquaternion of the 
form a + J7 with a, 7 £ R, is isomorphic to a complex number. Instead the complex conjugate must 
be seen as a fundamental operation. Note very carefully that a biquaternion does not commute with its 
complex conjugate, a surprising result that we examine in Proposition 1 below. 
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Of course, we can apply both conjugates [45] (we call this a total conjugate, but the term biconjugate 
has also been suggested in [12], and we consider it apposite). It is not difficult to show the following 
results: 

¥ = q* 

(pq) = p q 

In terms of the various representations above, the total or biconjugate is: 

q* = w r - x r i - y r j - z r k - I(w l - x t i - yj - z t k) =cfc- 1% = S(?)* ~ v (?)* ( 9 ) 
However, the ramifications of non-commutativity are deep, as shown by the following Proposition. 

Proposition 1. A biquaternion does not, in general, commute with its complex conjugate: reversing the 
order of the product yields a complex conjugate result. 

Proof. Represent an arbitary biquaternion in complex form I, as given in Table 2: q = q r + Iq^. Then 
the two products of q with its complex conjugate q* are: 

qq* = {qr + iqi) {q r - lq%) = qf + qf + i{q t q r - q r qi) 
q*q = (q r - Iqi) (qr + iqi) = qf + qf - i(q%Qt - q r q,) 

The results are a complex conjugate pair, and therefore differ. □ 
Note that important exceptions to Proposition 1 are: 

• Quaternions, that is biquatcrnions with q t = (trivial). 

• Imaginary biquaternions, with q r = 0. In this case the product is qf, regardless of order. 

• Biquatcrnions with real and imaginary parts that commute (co-planar, with a common axis). In 
this case the imaginary part of the product vanishes and the result is qf + qf, regardless of order. 

The Hamilton and complex conjugates correspond to negation of the argument of the exponential in 
the two polar forms in Table 2. In the 'Hamilton' polar form, negating the argument of the exponential 
negates the sine term, which corresponds to the vector part of q, and thus yields the classical Hamilton 
conjugate. In the 'complex' polar form, negating the argument of the exponential again negates the sine 
term, but in this case the sine term corresponds to the imaginary part of the quaternion, and thus it is 
the complex conjugate that is obtained. A third possibility would be a polar form that would correspond 
(in the sense just outlined) to the total conjugate. However, this would require the sine term to represent 
all the components in equation 9 except w r , which would be represented in the cosine term. This appears 
improbable. 

The three types of conjugate allow us to construct formulae for extracting the geometric components 
of a biquaternion, as defined in Table 3. The formulae in Table 4 are obtained by substitution from four 
formulae for the scalar /vector and real/imaginary parts of a quaternion based on sums and differences 
of the Hamilton or complex conjugates. Unlike the formulae in equation 2 however, which are based 
on involutions, and therefore require multiplication and addition only, these formulae contain complex 
conjugates, which are primitive operations that cannot be expressed in terms of multiplications and 
additions. Similarly, and fairly trivially, by combining the formulae in equation 2 with formulae based 
on sums and differences of complex conjugates, one can obtain formulae for extracting any of the eight 
real components of a biquaternion 9 . For example, that for the real part of the scalar part is w r — 
§(? + q* ~ iqi - iq*i - jqj - jq*j - kqk - kq*k). 

9 The utility of such formulae is not in numeric computation, where explicit access to the components is a much better 
approach, but in algebraic manipulation, where the formulae may allow an algebraic solution to an otherwise seemingly 
difficult derivation. 



8 



Scalar 


B(S(?)) 


1(9 + 






Bivector 


R(V(g)) 


1(9- 




-n 


Vector 


9f(V( 9 ))J 


1(9- 


9-9* 




Pscudoscalar 


9f(S(g))J 


l(<? + 


q-q* 


-r) 



Table 4: Formulae for the geometric components of a biquaternion. (See also Table 3.) 



3.2 Inner Product 

The definition of the inner product in the classic quaternion case was given by Porteous [32, Prop. 10.8, 
pl77] and we take this as the basis for defining the inner product of two biquaternions, since it is 
consistent with the quaternion case, and also yields the semi-norm (see the next section) in the case of the 
inner product of a biquaternion with itself. The same formula is given by Synge [43, Equation 3.8, p9] in 
the context of a discussion of basic properties of biquaternions. Ward, in his discussion of biquaternions 
and the Minkowski metric [45, §3.3, pll4], utilises the same definition, but he also discusses other 
definitions of the inner product (see [45, §3.2, pl09]). The definitions of Porteous and Synge are: 

(p,q) = 2 (pi + qp) = \ (pq + qp) (10) 

where the overbar represents a (Hamilton or quaternion) conjugate. The result of this expression will 
have zero vector part (since the vector parts of the two terms inside the parentheses cancel). In general, 
the scalar part of the result, and therefore the inner product, will be complex. The inner product may 
also be defined in terms of a simple clcmcntwisc product of the two quaternions, as can be seen by 
expanding out equation 10. If we let p = W p + X p i + Y p j + Z p k and q = W q + X q i + Y q j + Z q k, then 

{p,q) = W p W q + X p X q + Y p Y q + Z p Z q 

which is, of course, complex. (The result will be equal to the scalar part of the quaternion with zero 
vector part resulting from equation 10.) 

As already discussed in § 3, the use of the term inner product here does not imply that all the usual 
properties of an inner product will be satisfied. Ward [45, §3.3, pll5] states the following properties 
that are satisfied by the inner product defined in equation 10, and these are easily verified: 

(p,q) = (q,p) , p,qeM 

(p, q + r) = (p, q) + (p, r) , r G B 

a (p, q) = (ap, q) = (p, aq) , aeC 

Conventionally, but not always [31, See: Scalar product], the inner product is positive definite, that is 
greater than or equal to zero, or non-negative, which cannot apply to a complex-valued inner product. 

Classically, the inner product of two vectors is given by |i>i| |«2|cos0 where 9 is the angle between 
the two vectors. If extended to quaternions, it is not difficult to see that the angle is that between 
the two quaternions in 4-space, in the common plane defined by the two quaternions. When we extend 
this concept to biquaternions, the angle between them, and their moduli, must be complex, in general, 
so the geometric interpretation of the inner product is slightly more difficult. In the quaternion case, 
orthogonality arises from the angle between the two quaternions (having a vanishing cosine), but in the 
biquaternion case there is the additional possibility that one or both moduli are zero, resulting in a 
vanishing inner product. Expressing the inner product of two biquaternions in terms of inner products 
of the real and imaginary parts makes it possible to understand what the inner product represents 
geometrically. Representing p as follows: p — p r + Ipi and similarly for q, the inner product may be 
expanded as: 

(p,q) = (p r> q r ) - (Pi,q%) +i({p r ,Qi) + (Pi,qr) ) 

Certain special cases are apparent after careful inspection: 



9 



• (Pi <?) = indicates that the two biquaternions are orthogonal or 'perpendicular', and it requires 
the real and imaginary parts of the inner product to vanish separately. The detailed conditions 
required for this are many, since the four inner product components (two real, two imaginary) may 
have positive or negative real values and can cancel in several ways. We list here four different 
ways in which the scalar product can vanish: 

The strongest constraint is when all four real and imaginary parts of the two biquaternions are 
mutually orthogonal. A simple example is p = and q = j + Ik, but it is easy to construct 

a more general example starting from a random quaternion. Let p\ be a randomly chosen 
(pure or full) quaternion. Then let P2 = pii, P3 = pij and p^ = p\k. The four quaternions 
thus constructed have the same four numeric components permuted with sign changes in such 
a way that any two will have a vanishing inner product 10 . Two biquaternions p = p\ + Ipi 
and q = P3 + Ipi (or any other permutation) will have (p, q) = 0. This pair of biquaternions 
will be found to be divisors of zero (see § 3.5), but scaling the real and imaginary parts with 
different scale factors will result in biquaternions that are not divisors of zero, but are still 
orthogonal (orthogonality not being dependent on norm). Choosing the initial random value 
to be pure does not result in the pair of biquaternions being pure, because of permutation of 
the components. 

Weaker constraint I 

1. The real parts of the two biquaternions are orthogonal, and the imaginary parts of the two 
biquaternions arc orthogonal, hence the real part of the scalar product vanishes because 
it is the difference of two vanishing scalar products. 

2. The real part of each biquatcrnion is not orthogonal to the imaginary part of the other. 
However, the scalar products of the real part of one biquaternion with the imaginary part 
of the other have the same values, but opposite signs. This means the imaginary part of 
the scalar product of the two biquaternions vanishes by cancellation. 

A pair of biquaternions satisfying this constraint can be constructed from an arbitrary bi- 
quaternion p, by choosing q = pi (and two others also orthogonal to the first may be con- 
structed by choosing q = pj or q = pk) 11 . 

Weaker constraint II 

1. The real parts of the two biquaternions are not orthogonal, and neither are the imaginary 
parts, but the scalar product of the two real parts has the same value and sign as the 
scalar product of the two imaginary parts. This means the real part of the scalar product 
of the two biquaternions vanishes by cancellation. 

2. The real part of each biquaternion is orthogonal to the imaginary part of the other, hence 
the imaginary part of the scalar product of the two biquaternions vanishes because it is 
the sum of two vanishing scalar products. 

A pair of biquaternions satisfying this constraint can be constructed from an arbitrary bi- 
quatcrnion p, by choosing q = pii (and ditto, mutatis mutandis.). 

The weakest constraint is when none of the real and imaginary parts of the two biquaternions 
are orthogonal, but the real and imaginary parts of the scalar product of the two biquaternions 
vanish by cancellation. A pair of biquaternions satisfying this constraint can be constructed 
from an arbitrary biquatcrnion p, by choosing q = p(i + Ij) or similar. Note that this results 
in q being a divisor of zero (see § 3.5) because (i + Ij) is a divisor of zero. 

• (p,q) is real, that is with zero imaginary part. There is a trivial case: the two biquaternions 
may be imaginary with zero real parts. Otherwise, the imaginary part of the inner product can 
vanish in two ways, as in the weaker cases above of (p, q) = 0: either the imaginary part vanishes 
by cancellation, or it vanishes because the real part of each biquatcrnion is orthogonal to the 
imaginary part of the other. It is easy to construct a biquaternion to satisfy the second condition 
using the technique outlined above for the strongest constraint: construct two pairs of orthogonal 
quaternions, and then construct two biquaternions using the components of these pairs. 

10 Three other permutations can be obtained by multiplication on the right, or by division on either side by i, j and k. 
11 As in the stronger case, it is also possible to multiply by i on the left, etc. 
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• (p, q) is imaginary, that is with zero real part. The conditions required for this are very similar to 
the previous case, with the appropriate changes. A trivial case is where one biquaternion is real 
and the other is imaginary. 

Clearly, from the above analysis, orthogonality of biquaternions is not as simple as orthogonality of real 
quaternions where a geometrical interpretation in terms of a plane and a real angle is straightforward. 
In the biquaternion case, interpreting the inner product (p, q) as \p\ \q\ cos O, where the two moduli and 
the angle are complex, is not at all obvious. This remains a topic for further work. 

3.3 Semi-norm 

It is possible to define more than one norm for the biquaternions (sec for example [45, §§3.2 and 3.3]). 
In this section we consider the generalization of the quaternion norm \\q\\ ~ w 2 + x 2 + y 2 + z 2 , 1 1 g 1 1 € R to 
the biquaternion case by allowing the four Cartesian components to become complex. Conventionally a 
norm is real and positive definite, that is non-negative, but in the case of biquaternions, this convention 
must be relaxed because the norm has a complex value. Additionally, the existence of divisors of zero 
(see § 3.5) means that a non-zero biquaternion can have a vanishing norm. For this reason, we use the 
term semi-norm. A semi-norm is a generalization of the concept of a norm, with no requirement that 
the norm be zero only at the origin (of a vector space) [3, See: semi-norm]. As with our use of the term 
inner product in § 3.2 for a complex-valued generalization of the inner product, we are here extending 
the term semi-norm to a complex-valued quantity analogous to a norm, with the additional property 
of vanishing in the case of divisors of zero. We also use the term modulus for the square root of the 
semi-norm - it is also of course complex, in general. 

The semi-norm can be real and negative in special cases, as well as purely imaginary. The semi-norm 
can be defined in terms of the inner product of a biquaternion with itself, that is |jg|| = (q,q) or directly 
in terms of the four complex components: ||g|| = \q\ = W 2 + X 2 + Y 2 + Z 2 . Although the semi-norm 
is complex-valued the result given by Coxeter for quaternions [6, §2] still holds: \\q\\ = \q\" = qq = qq. 
This is, of course, a special case of the formula for the inner product, in this case of q with itself (equation 
10). 

Lemma 1. [35, 36, Lemma 1] Let q = q r + Iqi be a non-zero biquaternion with real part q r £ H and 
imaginary part q,; € H. The real part of \\q\\ is equal to the difference between the norms of q r and qi, 
and the imaginary part of \\q\\ is equal to twice the inner product of q r and q^. 

Proof. We express the semi- norm of q as the product of q with its quaternion conjugate: ||q|| = qq. 
Writing this explicitly: 

IMI = (<?r + Iqi)(q r + Iqi) = (q r + iQijiW + iQi) 

and multiplying out we get: 

1 1 q 1 1 = q r q~r - liQi + Hlrql. + qiCfr) 

The real part of H^H can be recognised as — \\q.i\\ from Coxcter's result [ , § 2]. The imaginary part 
of \\q\\ is twice the inner product of q r and qi as given by Porteous [32, Prop. 10.8, pl77] and equation 
10. □ 

There are some special cases of the result given in Lemma 1: 

• biquaternions with perpendicular real and imaginary parts have a real semi-norm (because the 
imaginary part vanishes), but the semi- norm may be negative if the norm of the imaginary part 
exceeds the norm of the real part; 

• biquaternions with real and imaginary parts with equal norms have an imaginary semi-norm (be- 
cause the real part vanishes). 

• imaginary biquaternions (with q r = 0) have a negative real semi-norm, and therefore an imaginary 
modulus. 
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• biquaternions with perpendicular real and imaginary parts with equal norms have a vanishing 
semi-norm. This important special case is covered in more detail in § 3.5. 

The semi-norm is invariant under quaternion conjugation, that is \\q\\ = \\q\\, but not under complex 
conjugation: \\p\\ = \\p*\\ ■ This can be seen easily from Lemma 1, since quaternion conjugation does not 
affect the norms of the real and imaginary parts, nor their inner product. However, although complex 
conjugation docs not alter the norms of the real and imaginary parts, it does negate their inner product 
and therefore negates the imaginary part of the semi-norm. 

The semi-norm as defined here (and by Ward [45, §3.3, pll5]) obeys the 'rule of the norms' 12 so 
that for any two biquaternions, the semi-norm of their product equals the product of their semi-norms: 
\\P c l\\ = WpW \\q\\- This is true even in the case where one or both of p and q is a divisor of zero, with 
vanishing semi-norm, as discussed in § 3.5 (but the result will be zero, computed as the semi-norm of 
the product, or the product of the semi-norms). Hence B is a normed algebra using the definition of the 
semi-norm used in this section. The non-zero biquaternions that have zero semi-norm (the divisors of 
zero) lack a multiplicative inverse. 

A norm is usually required to satisfy the following three properties [3, See: norm] (caution 13 ): 

IM| = |A||MI (ii) 

1!^ + y\\ < ||x|| + \\y\\ (Triangle inequality) 

The first of these properties is easily seen to be satisfied for all q € B because negating a complex number 
does not change its square, and therefore does not alter the sum of the squares of the four components 
of a biquaternion, nor the square root of the sum of the squares. The second property requires rather 
more careful analysis, but when the notation is correctly interpreted, it can be seen to be satisfied by the 
square root of the semi-norm of a biquaternion. The quantity A is usually taken to be a scalar, that is in 
the biquaternion context, a complex number. If we write the expression explicitly in biquaternion terms, 
using the Cartesian form, we have: |Ag| = ^(XW) 2 + (XX) 2 + (XY) 2 + (XZ) 2 = y/j?\q\. Taking the 
square root of the square of a scalar quantity conventionally yields the positive square root, or absolute 
value. In the case of a complex A, taking the square root of the square yields a result in the right 
half-plane, rather than the absolute value. Thus to decide whether the biquaternion semi-norm satisfies 
the property as conventionally stated requires some deeper understanding of the reason for the property. 
The triangle inequality, as has already been noted in § 3 is not applicable to the biquaternion semi-norm, 
because ordering is not defined for complex numbers. 

Finally we note that it is possible to define a real norm and modulus for the biquaternions. Giirlcbcck 
and Sprossig [13, Lemma 1.30] state that there is a unique real norm satisfying \\pq\\ = \\p\\ \\q\\ and 
IIm^oII = ImI ;A* G R where Co is a basis element (Pauli matrix). This norm is simply the absolute value 
of the square root of the semi-norm, in other words, r in equations 4 and 5, or \Q\ in Table 2. 

3.4 Roots of -1 

This section draws on and builds on results first presented by Sangwine in [33, 34]. These results have 
been generalized to Clifford algebras by Hitzer and Ablamovicz [25, 26] where the concept is aptly referred 
to as geometric roots of —1. 

Biquaternion roots of -1 appear as £ in Table 2, in the 'complex' and 'Hamilton polar' forms. 

Theorem 1. Any non-zero quaternion or biquaternion with a non-vanishing modulus, divided by its 
modulus has a unit real norm. 

Proof. Let q be an arbitrary quaternion or biquaternion and let p = q/ \q\. For any quaternion or 
biquaternion x we have ||x|| = xx and |x| = ^/jjx][. Therefore p = q/y/qlj and we can write the norm of 

12 The 'rule of the norms' holds also for the modulus (the square root of the norm), which is what the term 'norm' refers 
to in many sources in the literature. 

13 The properties given in equation 11 are stated using the notation used in the cited reference which does not conform 
to the usage elsewhere in this paper. In particular, in equation 11 means what we call the modulus or square root of 
the norm in this paper. 
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p as: 



'1 



Q 



\\P\\ =PP 



Since y/qlj is complex, it is unaffected by the quaternion conjugate and we can simplify this to: 
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? = q_q 

qq^/qq qq 



= 1 



□ 



Theorem 2. Any pure quaternion or biquaternion with a non-vanishing modulus, divided by its modulus, 
is a root of —I. 

Proof. Let q be an arbitrary quaternion or biquaternion with \q\ ^ 0. Then (q/ \q\) 2 = q 2 /qq = q/q. If 
we restrict q to have zero scalar part, the conjugate reduces to negation and we have (q/ \q\) = qj — q = 
-1. □ 

Theorem 2 is simply a restatement of the well-known fact that unit pure quaternions are roots of 
— 1 [17, § 167, pl79], but in the case of biquaternions, the ramifications of the result are not so obvious. 
Theorem 2 is also easily demonstrated using the Cartesian form: 



Xi + Yj + Zk 

VX 2 + Y 2 +1* 



-X 2 + XY{ij + ji) + XZ(ik + ki) -Y 2 + YZ(jk + kj) - Z 2 



X 2 + Y 2 + Z 2 



- X 2 + Y 2 



Z 2 



X 2 +Y 2 + Z 2 



(12) 



It was shown in [34, Theorem 2.1] 14 that any pure biquaternion £ satisfying the constraints _L 
and ||9?(£)|| — 11^(011 = 1 is a root of —1 and that no other biquaternions are roots of —1. It follows 
therefore that dividing an arbitrary pure biquaternion by its (complex) modulus produces a result that 
satisfies the constraints stated. Since these constraints are by no means obvious from equation 12, 
it is interesting to verify them directly. To verify the constraints stated, we take an arbitrary pure 
biquaternion p, divide it by its (complex) modulus and show that the result satisfies the constraints 
stated. Dividing p by its modulus gives a root of —1: £ = pj \p\. We have to show that the real and 
imaginary parts of £ arc orthogonal, that is (5i(£), $i(£)) = Oj an d also that the difference between the 
norms of the real and imaginary parts of £ is 1. We can express the real and imaginary parts of £ as 
follows, by taking the sum and difference of £ with its complex conjugate: 



&(0 = o(£ + O> 3(0 



21 



(€-€*) 



The inner product of two quaternions is given by equation 10, but since we are dealing with pure 
quaternions, the conjugates reduce to negation, and the inner product simplifies to: 



1 



(w, v) 



{uv + vu) 



(13) 



Substituting the results above for the real and imaginary parts of we obtain the following: 



= -gj [(€ + **) (€-€*) + (€ -€*)(* + €*)] 



l 

"81 



*Also in [33, Theorem 1]. 
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Since £ is a root of —1 by Theorem 2, its complex conjugate must also be a root of minus one 15 , thus 
i 2 = C 2 = -1 and therefore = as expected. 

To show that the real and imaginary parts of £ have norms differing by 1, we can use a similar 
approach. Starting with the fact that the norm is the inner product of a quaternion with itself, as stated 



in § 3.3, the formula in equation 13 simplifies to 



when u is pure and v = u. Thus the difference 



between the norms of the real and imaginary parts of £ is given by: 



(OIM|s(OII = -^) 2 + 3(£) 2 



21 



(€-€*) 



e + e - «* - en 



which is 1 because £ and its conjugate are roots of — 1. 



3.5 Divisors of zero 

It is possible for the semi-norm, \\q\\, of a biquatcrnion to be zero, even though q ^ 0. In modern 
terminology, the biquatcrnions with vanishing semi-norm are better known as divisors of zero. The 'rule 
of the norms' given in § 3.3 shows that multiplying an arbitrary biquatcrnion by a divisor of zero (with 
vanishing semi-norm) yields a result which is also a divisor of zero (with vanishing semi- norm) . 

The conditions for the semi-norm to vanish were discovered by Hamilton [17, Lecture VII, §672, 
p669] 16 — if we represent a non-zero biquatcrnion in the form: q = q r + Iqi then its norm is zero 
iff \\q r \\ = \\qi\\ and (q r ,qi) = 0, that is, the real and imaginary parts must have equal norms and be 
perpendicular in 4-space. This result follows easily from Lemma 1 in § 3.3. The conditions stated are 
equivalent to the condition W 2 + X 2 + Y 2 + Z 2 = 0, which can be expanded into a pair of simultaneous 
equations each equivalent to one of the two conditions above, by equating the real and imaginary parts 
to zero separately [35, 36, Proposition 1]. Notice that the conditions required for a vanishing semi- 
norm may be satisfied in the case of a pure biquatcrnion (the real and imaginary components are then 
perpendicular in 3-spacc). It follows from these conditions that any non-zero quaternion or any non-zero 
imaginary biquaternion must have a non-zero norm, since the constraint of real and imaginary parts 
having equal norm is violated. 

Sangwine and Alfsmann have shown in [35, 36] that every biquaternion divisor of zero is one of: 

• an idempotent as defined in the next section in Theorem 3; 

• a complex multiple of an idempotent (dividing the divisor of zero by twice its scalar part yields an 
idempotent); 

• a nilpotent as defined in §3.5.2. 



3.5.1 Idempotents 

An idempotent is a value that squares to give itself. An example of a biquaternion idempotent is given 
in [12] (quoting a paper of Lanczos) but no general case is given. We present here a derivation of the set 
of idempotents in B, drawing on results first presented in [35, 36] by Sangwine and Alfsmann. 

Theorem 3. [35, 36, Theorem 2] Any biquaternion of the form q = | ± g£7 , where £ £ B is a root of 
— 1, is an idempotent. There are no other idempotents in B. 

Proof. The proof is by construction from an arbitrary biquatcrnion represented in the form q = A + £B 
{A,B € C). Squaring and equating to q gives: 

q 2 = A 2 - B 2 + 2AB£ = A + £B =q 

15 This is easily shown using the 'complex' form £ = £ r + since the imaginary part must be zero in both £ 2 and £* 2 
if they are each equal to —1. 

16 In [17, Lecture VII, §673, p671], Hamilton rather quaintly refers to an evanescent tensor — a tensor being what we 
here call a semi-norm. 
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Equating coefficients of 1 and £ (that is, equating the reai and imaginary parts with respect to £) gives: 



A = A 2 - B 2 



B = 2AB 




The roots of —1 in B include two trivial cases, as noted in [34, Theorem 2.1]. These are ±7 which 
yields the trivial idempotents q = and q = 1; and ±/x where \x is a real pure quaternion root of —1. 

Since every biquaternion idempotent, q, is a solution of q 2 = 0, it must also be a solution of q(q~l) = 
0. Therefore every biquaternion idempotent is also a divisor of zero and must satisfy the conditions stated 
in §3.5. This is demonstrated in detail in [35, 3G, §4, Proposition 2]. 

All non-pure biquaternion divisors of zero have been shown to be of the form p = aq where a is a 
non-zero complex number, and q is a biquaternion idempotent [35, 36, § 4, Theorem 3]. Further, dividing 
a non-pure biquaternion divisor of zero by twice its scalar part yields an idempotent [35, 36, § 4, Corollary 



3.5.2 Nilpotents 

A nilpotcnt is a quantity whose square vanishes. Hamilton discovered nilpotcnt biquatcrnions and was 
aware that all such biquatcrnions were divisors of zero [17, Lecture VII, §674, pp671-3]. 

We first present the case of a pure biquaternion with vanishing semi-norm and show a result due to 
Hamilton [17, Lecture VII, § 672, p669] that the square of such a biquaternion vanishes. 

Theorem 4. Let p = p r + Ipt be a non-zero pure biquaternion with vanishing semi-norm. Then p 2 = 0. 

Proof. From the fact that p has a vanishing semi-norm, ||p r || = ||pj||, therefore we may divide by their 
common norm, and obtain pj ||p r || = \x + Ii> where (x and v are perpendicular unit pure quaternions. 



which vanishes because the squares of the two unit pure quaternions are —1 and the imaginary part 
vanishes because the products of two perpendicular unit pure quaternions changes sign when the order 



Sangwine and Alfsmann, in [35, 36, § 5, Lemma 1] have shown that all nilpotent biquatcrnions are 
pure, and in [35, 36, § 5, Proposition 3] that all nilpotent biquatcrnions are divisors of zero. It follows 
that all nilpotent biquaternions have real and imaginary parts with equal norm and therefore that any 
nilpotcnt biquaternion can be normalised to the form //, + Iv as in Theorem 4. We see therefore that all 
nilpotcnt biquaternions are constructed from a pair of mutually orthogonal unit pure quaternions. 

A nilpotent biquaternion can be written in the form £B where £ is a root of —1 and B is complex 
(the square root of the semi- norm: B = V 1 X 2 + Y 2 + Z 2 ). Since a nilpotcnt biquaternion is a divisor of 
zero, its semi-norm is zero. This means that B is zero, even though £B is not. A consequence is that it 
is not possible to compute the 'axis' of a nilpotcnt biquaternion by dividing by B. 

3.6 The biquaternions as a geometric algebra 

In this section we discuss equivalence of the components of a biquaternion to the elements of a geometric 
algebra [24, 40, 23, 30, 42] or Clifford algebra. Jaap Suter's tutorial paper [42] is particularly recom- 
mended. A full discussion of the biquaternions as a geometric algebra is a large and non-trivial topic, 
and will be left to a later paper or papers. 

Hamilton himself, in his Lectures on Quaternions [17, Lectures I, II, III] was interested in an algebra 
of points, lines and planes, and clearly was motivated by a desire to have an algebra applicable to 3- 
dimcnsional geometry The modern notion of a geometric algebra uses different language, and some of 



1]. 



Then: 




of the product is reversed. 



□ 
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the concepts are more sophisticated, but essentially the aim is the same: to represent and manipulate 
geometric objects algebraically. 

A geometric algebra contains four types of element, as already introduced in §3 and Tables 3 and 4: 
scalars, vectors, bivectors and pseudoscalars. Scalars are quantities without geometric form (or points), 
vectors represent directed magnitudes (or lines), bivectors represent directed areas, and pseudoscalars 
represent signed volumes. The product of two perpendicular vectors is a bivector, representing the area 
swept out by one vector when moved along the other [42, §2.1]. This is why vectors arc identified 
with pure imaginary biquaternions. Given two perpendicular pure quaternions, fi and u, such that 
(fi, v) = 0, the two vectors Ifi and Iu have a product — fiv which is perpendicular to both fi and v. 
The negative sign is a consequence of squaring I. Changing the order of the two vectors changes the 
sign of the resulting bivector (it is oppositely directed), so the minus sign is of little consequence. It is 
usual in geometric algebra to have to choose the sense by convention. Here the rules of biquaternion 
multiplication define the convention for us. The product of three orthogonal vectors defines a volume, 
represented by a pseudoscalar. For example, (Ii)(Ij)(Ik) = I, as given by Ward [45, pll3]. 

There has been much confusion caused by the use of the term vector when in fact the objects under 
consideration were bivectors. In physics the concepts of axial and polar vectors have been used, causing 
further confusion, because the two terms suggest different types of vector, rather than fundamentally 
different types of quantity. Even the present authors, in their earlier works, used the term vector to refer 
to quantities which are now clearly seen to be bivectors. Jaap Suter in his Geometric Algebra Primer 
states clearly: 

A quaternion is a scalar plus a bivector. 

and he is right [42, p50]. Unfortunately, the scalar/vector part terminology commonly applied to quater- 
nions has this wrong, and we are stuck with it. 

In geometric algebra a composite quantity consisting of two or more of the four types of quantity 
(scalar, vector, bivector, pseudoscalar) is called a multivector. A biquaternion may contain all four 
types of geometric element, and except in degenerate cases (such as a pure imaginary biquaternion), 
corresponds to the concept of a multivector. 

Once we start to multiply multivectors, the geometric interpretation becomes more difficult, although 
the rules of multiplication are well-understood and easily described: the full biquaternion product can 
be expressed, exactly as in the quaternion case, as: 

pq = S(p) S(q) + S(p) V(q) + S(q) V(p) + V(p) V(q) 

and the last term can be written as V(p) ■ V(g) + V(p) x V(q) or (V(p) , V(g)) + V(p) x V(q), that is 
the sum of a dot or scalar or inner product, and a cross product. A deeper analysis of the biquaternions 
as a geometric algebra requires further work. 

Although multivectors in geometric algebra are regarded as composites of scalar, vector, bivector and 
pseudoscalar, a biquaternion is not constructed in this way. Instead it has four complex components, 
but as we have seen in § 3, it is possible to interpret these in multiple ways, for example as a complex 
number with quaternion components. It is also possible to consider the biquaternions as 8-dimcnsional 
quantities with basis (l,il,jl,kl,i,j,k,l) where we place the vector basis first. The multiplication 
table is then easily derived from the rules of quaternion multiplication (Table 5). 

From the rules of quaternion multiplication, it is fairly easy to draw up a table showing how the 
geometric components of a biquaternion multiply together. Table 6 shows the result. The four centre 
entries apply in the general case, but in the specific cases where the bivector /vector are perpendicular, 
the scalar or pseudoscalar part of the product will be zero. Multiplication by the pseudoscalar is known 
as the dual operation: it maps a bivector into a vector, and a scalar into a pseudoscalar or vice versa. 

3.6.1 The wedge or outer product 

In geometric algebra, the algebraic or geometric product of two vectors is the sum of the scalar product 
and the wedge (A) or outer product (also known as the cross product in three dimensions): 

pq = p ■ q + p Aq 
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Tabic 5: Biquatcrnion basis multiplication tabic. 
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Table 6: Multiplication table for components of a biquaternion multivector 

It is possible to define a wedge product for quaternions and biquatcrnions using a formula widely used 
in geometric algebra [42, Equation 3.10]: 

PA q = - (pq - qp) 

This gives the cross product of the vector parts (because the components of the products involving 
the scalar parts, and the dot product of the vector parts, commute and therefore cancel). If applied 
to a pair of vectors (i.e. pure biquaternions with zero real part), the result is a bivector (i.e. a pure 
biquaternion with zero imaginary part). Applied to a pair of bivectors (pure quaternions) the result is 
another bivector. Applied to a vector and a bivector, the result is a vector. 

4 Further work 

There is much further interpretation needed of the biquaternions as a geometric algebra. At present 
formulae are lacking for vector/bivector products to yield pseudoscalars and for generalisations of the 
wedge product. 

The interpretation of complex angles in the polar form and in the context of the inner product requires 
further work. 
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